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T-H We survey an area of recent development, relating dynamics to theo- 

retical computer science. We discuss the theoretical limits of simulation 

' _ ' and computation of interesting quantities in dynamical systems. We will 

focus on central objects of the theory of dynamics, as invariant measures 

(— H and invariant sets, showing that even if they can be computed with arbi- 

(— j trary precision in many interesting cases, there exists some cases in which 

4-^ they can not. We also explain how it is possible to compute the speed 

a of convergence of ergodic averages (when the system is known exactly) 
and how this entails the computation of arbitrarily good approximations 
' ' of points of the space having typical statistical behaviour (a sort of con- 

structive version of the pointwise ergodic theorem). 
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1 Introduction 

The advent of automated computation led to a series of great successes and 
achievements in the study of dynamical problems. 

The use of computers and simulations allowed to compute and forecast the 
behavior of many important natural phenomena and, on the other hand, led to 
the discovery of important general aspects of dynamics. 

This motivates the huge work that was made by hundred of scientists to 
improve "practical" simulation and computation techniques. 

It also motivates the study of the theoretical limits of simulation and com- 
putation techniques, and the theoretical understanding of related problems. 

In this paper we want to focus on some of these aspects related to rigorous 
computation and simulation of (discrete time) dynamical systems. 

The simulation and investigation of dynamics started with what we call 
the "naive" approach, in which the user just implements the dynamics without 
taking rigorous account of numerical errors and roundings. Then he "looks" to 
the screen to see what happens. 

Of course, the sensitivity to initial conditions, and the typical instability of 
many interesting systems (to perturbations on the map generating the dynam- 
ics) implies that what it is seen on the screen could be completely unrelated to 
what was meant to be simulated. 

In spite of this, the naive approach turns out to work "unreasonably" well in 
many situations and it is still the most used one in simulations. The theoretical 
reasons why this method works and its limits, in our opinion are still to be 
understood (some aspects have been investigated in [36], [41], [28] [7, 8] e.g.). 

On the opposite side from the naive approach, there is the "absolutely rig- 
orous" approach, which will be the main theme of this paper: the user looks 
for an algorithm which can give a description of the object which is meant to 
be computed, up to any desired precision. 

In this point of view, for example the constant e is a computable number 
because there is an algorithm that is able to produce a rational approximation 
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of e at any given precision (for example finding the right m and calculating 
^™ ^ such that the error is smaller than requested). 

In an "absolutely rigorous" simulation (computation) the initial point or 
some initial distribution is supposed to be known up to any approximation and 
the transition map is also supposed to be known up to any accuracy (suitable 
precise definitions will be given below) . This allows the evolution of the system 
to be simulated with any given accuracy, and the question arise, if interest- 
ing and important objects related to the dynamics (invariant sets or invariant 
measures e.g.) can be computed from the description of the system (or maybe 
adding some additional information). 

In this paper we would like to give a survey of a group of results related 
to these computational aspects, restating and updating them with some new 
information. We will see that in many cases the interesting objects can be 
computed, but there are some subtleties, and cases where the interesting object 
cannot be computed from the description of the system or cannot be computed 
at all (again, up to any given precision). In particular, this happen case for the 
computation of invariant measures and Julia sets. 

Hence, these results set theoretical limits to such computations. 

Computing invariant measures. 

An important fact motivating the study of the statistical properties of dy- 
namical systems is that the pointwise long time prediction of a chaotic system is 
not possible, whereas, in many cases, the estimation or forecasting of averages 
and other long time statistical properties is. This often corresponds in math- 
ematical terms to computing invariant measures, or estimating some of their 
properties, as measures contain information on the statistical behavior of the 
system {X, T) and on the potential behavior of averages of observables along 
typical trajectories of the system (see Section 3). 

An invariant measure is a Borel probability measure fi on X such that for 
each measurable set A it holds n{A) — ^{T~^(A)). They represent equilibrium 
states, in the sense that probabilities of events do not change in time. 

Rigorously, compute an invariant measure means to find an algorithm which 
is able to output a description of the measure (for example an approximation 
of the measure made by a combination of delta measures placed on "rational" 
points) up to any prescribed precision. 

We remark that once an interesting invariant measure is computed, it is pos- 
sible to deduce from it several other important information about the dynamics: 
Lyapunov exponents, entropy, escape rates, etc... For example, in dimension 
one, once an ergodic invariant measure fi has been approximated, the Lyapunov 
exponent can be estimated using the formula = log2 T'dfi, where T' de- 
notes the derivative of the map generating the dynamics. In higher dimensions, 
similar techniques can be applied (see e.g. [21] for more examples of derivation 
of dynamical quantities from the computation of the invariant measure) . 

Before giving more details about the computation of invariant measures, we 
remark that, since there are only countably many "algorithms" (computer pro- 
grams), whatever we mean by "approximating a measure by an algorithm" would 
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imply that only countable many measures can be computed whereas, in general, 
a dynamical system may have uncountably many invariant measures (usually an 
infinite dimensional set). So, a priori most of them will not be algorithmically 
describable. This is not a problem because we should put our attention on the 
most "meaningful" ones. An important part of the theory of dynamical systems 
is indeed devoted to the understanding of "physically" relevant invariant mea- 
sures. Informally speaking, these are measures which represent the asymptotic 
statistical behavior of "many" (positive Lebesgue measure) initial conditions 
(see Section 3 for more details). 

The existence and uniqueness of physical measures is a widely studied prob- 
lem (see [47]), which has been solved for some important classes of dynamical 
systems. These measures are some of the good candidates to be computed. 

The more or less rigorous computation of such measures is the main goal 
of a part of the literature related to computation and dynamics. A main 
role here is played by the transfer operator induced by the dynamics. Indeed, 
the map T defining the dynamics, induces a dynamics Lt on the space of 
probability measures on X, Lt ■ PM{X) — )■ PM{X). Lt is called the transfer 
operator associated to T (definition and basic results about this are recalled 
in Section 3). Invariant measures are fixed points of this operator. The main 
part of the methods which are used to compute invariant measures deals with 
suitable finite dimensional approximation of this operator. In Section 3 we will 
review briefly some of these methods, and some references. We then consider the 
problem from an abstract point of view, and give some general result on rigorous 
computability of the physical invariant measure. In particular we will see that 
the transfer operator is computable up to any approximation in a general context 
(see Thm 15). The invariant measure is computable, provided we are able to 
give a description of a space of "regular" measures where the physical invariant 
measure is the only invariant one (see Thm 16 and following corollaries). 

We will also show how the description of the space can be obtained from the 
one of the system in a class of examples (piecewise expanding maps) by using 
the Lasota-Yorke inequality. 

After such general statements one could conjecture that all computable dy- 
namical systems should always have a computable invariant measure. We will 
see that, perhaps surprisingly, this is not true. Not all systems that can be de- 
scribed explicitly (the dynamics can be computed up to any prescribed approxi- 
mation) have computable invariant measures (see Section 3.5). The existence of 
such examples reveals some subtleties in the computation of invariant measures. 

To further motivate these results, we finally remark that from a technical 
point of view, computability of the considered invariant measure is a requirement 
in several results about relations between computation, probability, randomness 
and pseudo-randomness (see Section 6 and e.g. [3, 23, 24, 25]) 

Computability in Complex Dynamics 

Polynomial Julia sets have emerged as the most studied examples of fractal 
sets generated by a dynamical system. One of the reasons for their popularity 
is the beauty of the computer-generated images of such sets. The algorithms 
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used to draw these pictures vary; the naive approach in this case works by 
iterating the center of a pixel to determine if it Hes in the Juha set. There 
exists also more sophisticated algorithms (using classical complex analysis, see 
[42]) which work quite well for many examples, but it is well known that in 
some particular cases computation time will grow very rapidly with increase of 
the resolution. Moreover, there are examples, even in the family of quadratic 
polynomials, where no satisfactory pictures of the Julia set exist. 

In the rigorous approach, a set is computable if, roughly speaking, its image 
can be generated by a computer with an arbitrary precision. Under this notion 
of computability, the question arise if Julia sets are always computable. In a 
series of papers ([5, 4, 11, 12]) it was shown that even though in many cases 
(hyperbolic, parabolic) the Julia set is indeed computable, there exists quadratic 
polynomials which are computable (again, in the sense that all the trajectories 
can be approximated by an algorithm at any desired accuracy), and yet the 
Julia set is not. 

So we can not simulate the set of limits points on which the chaotic dy- 
namics takes place, but, what about the statistical distribution?. In fact, it was 
shown by Brolin and Luybich that there exists a unique invariant measure which 
maximizes entropy, and that this measure is supported on the Julia set. The 
question of whether this measure can be computed has been recently solved in 
[6], where it is proved that the Brolin-Lyubich measure is always computable. 
So that even if we can not visualize the Julia set as a spatial object, we can 
approximate its distribution at any finite precision. 

Computing the speed of convergence and pseudorandom points. 

In several questions in ergodic theory The knowledge of the speed of conver- 
gence to ergodic behavior is important to deduce other practical consequences. 
In the computational framework, the question turn out to be the effective 
estimation of the speed of convergence in the ergodic theorems^. From the 
numerical-practical point of view this has been done in some classes of systems, 
having a spectral gap for example. In this case a suitable approximation of the 
transfer operator allows to compute the rate of decay of correlations [22] [37] and 
from this, other rates of convergence can be easily deduced. 

Other classes of systems could be treated joining the above spectral ap- 
proach, with combinatorial constructions (towers, see [38] e.g.), but the general 
case need a different approach. 

In [2] it was shown that much more in general, if the system can be described 
effectively, then the rate of convergence in the pointwise ergodic theorem can be 
effectively estimated. We give in section 5 a very short proof of a statement of 
this kind (see Theorem 39 ) for ergodic systems, and show some consequences. 
Among these, a constructive version of pointwise ergodic theorem. If the system 
is computable (in some wide sense that will be described) then, it is possible to 
compute points having typical statistical behavior. Such points could be hence 
called pseudorandom points in the system (see Section 6). 

^Find a TV such that ^ foT" differs from J f dfi loss than a given error for each n> N. 
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Since the computer can only handle computable initial conditions, any sim- 
ulation can start only from these points. Pseudorandom initial conditions are 
hence in principle good points where to start a simulation. 

We remark that it is widely believed that naive computer simulations very 
often produce correct statistical behavior. The evidence is mostly heuristic. 
Most arguments are based on the various "shadowing" results (see e.g. [33] 
chapter 18). In this kind of approach (different from ours), it is possible to 
prove that in a suitable system every pseudo-trajectory, as the ones which are 
obtained in simulations with some computation error, is close to a real trajectory 
of the system. However, even if we know that what we see in a simulation is near 
to some real trajectory, we do not know if this real trajectory is typical in some 
sense. A limit of this approach is that shadowing results hold only in particular 
systems, having some strong hyperbolicity, while many physically interesting 
systems are not like this. In our approach we consider real trajectories instead 
of "pseudo" ones and we ask if there is some computable point which behaves 
as a typical point for the dynamics. 

Acknowledgement 1 We would like to thank The Abdus Salam International 
Centre for Theoretical Physics ( Trieste, IT) for support and hospitality during 
this research. 

2 Computability on metric spaces 

To have formal results and precise assumptions on the computability (up to any 
given error) of continuous objects, we have to introduce some concepts. 

We have to introduce some recursive version of open and compact sets, and 
characterize the functions which are well suited to operate with those sets (com- 
putable functions) . We explain this theory in this section trying the explanation 
to be to be as much as possible simple and self contained. 

2.1 Computability 

The starting point of recursion theory was to give a mathematical definition 
making precise the intuitive notions of algorithmic or effective procedure on 
symbolic objects. Several different formalizations have been independently pro- 
posed (by Church, Kleene, Turing, Post, Markov...) in the 30's, and have proved 
to be equivalent: they compute the same functions from N to N. This class of 
functions is now called the class of recursive functions. As an algorithm is al- 
lowed to run forever on an input, these functions may be partial, i.e. not defined 
everywhere. The domain of a recursive function is the set of inputs on which 
the algorithm eventually halts. A recursive function whose domain is N is said 
to be total. 

We now recall an important concept from recursion theory. A set _E C N 
is said to be recursively enumerable (r.e.) if there is a (partial or total) 
recursive function (p : N — > N enumerating E, that is E = {f{n) : n g N}. If 
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E ^ 9, if can be effectively converted into a total recursive function ip which 
enumerates the same set E. 

2.2 Algorithms and uniform algorithms 

Strictly speaking, recursive functions only work on natural numbers, but this 
can be extended to the objects (thought as "finite" objects) of any countable 
set, once a numbering of its elements has been chosen. We will use the word 
algorithm instead of recursive function when the inputs or outputs are inter- 
preted as finite objects. The operative power of algorithms on the objects of 
such a numbered set obviously depends on what can be effectively recovered 
from their numbers. 

More precisely, let X and Y be numbered sets such that the numbering of 
X is injective (it is then a bijection between N and X). Then any recursive 
function (/3 : N — > N induces an algorithm A : X Y. The particular case 
X = N will be much used. 

For instance, the set Q of rational numbers can be injectively numbered 
Q = {^Oi^ii---} in effective way: the number z of a rational a/b can be 
computed from a and 6, and vice versa. We fix such a numbering: from now 
and beyond the rational number with number i will be denoted by qi. 

Now, let us consider computability notions on the set M of real numbers, 
introduced by Turing in [44]. 

Definition 2 Let x be a real number. We say that: 

9 X is lower semi-computable if the set {i £ N : Qi < x} is r.e. 

• X is upper semi-computable if the set {i E N : Qi > x} is r.e. 

• X is computable if it is lower and upper semi- computable. 

Equivalently, a real number is computable if and only if there exists an algo- 
rithmic enumeration of a sequence of rational numbers converging exponentially 
fast to X. That is: 

Proposition 3 A real number is computable if there is an algorithm A : N ^ 
Q such that \A{n) - x\ < 2"" for all n. 

Uniformity. Algorithms can be used to define computability notions on many 
classes of mathematical objects. The precise definitions will be particular to 
each class of objects, but they will always follow the following scheme: 

An object O is computable if there is an algorithm 
A-.X-^Y 
which computes O in some way. 
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Each computability notion comes with a uniform version. Let (Oi)ieN be a 
sequence of computable objects: 

Oi is computable uniformly in i if there is an algorithm 
A-.NxX 

such that for all i, Ai := A{i, ■) : X computes Oj. 

For instance, the elements of a sequence of real numbers (.Xi)igN are uniformly 
computable if there is a algorithm ^ : N x N — )■ Q such that \A{i, n) — Xi\ < 2^" 
for all i, n. 

In other words a set of objects is computable uniformly with respect to 
some index if they can be computed with a "general" algorithm starting from 
the value of the index. 

In each particular case, the computability notion may take a particular name: 
computable, recursive, effective, r.e., etc. so the term "computable" used above 
shall be replaced. 

2.3 Computable metric spaces 

A computable metric space is a metric space with an additional structure allow- 
ing to interpret input and output of algorithms as points of the metric space. 
This is done in the following way: there is a dense subset (called ideal points) 
such that each point of the set is identified with a natural number. The choice of 
this set is compatible with the metric, in the sense that the distance between two 
such points is computable up to any precision by an algorithm getting the names 
of the points as input. Using these simple assumptions many constructions on 
metric spaces can be implemented by algorithms. 

Definition 4 A computable metric space (CMS) is a triple X = [X,d,S), 
where 

(i) {X,d) is a separable metric space. 

(ii) S = {sijigN is a dense, numbered, subset of X called the set of ideal 
points. 

(Hi) The distances between ideal points d{si,Sj) are all computable, uniformly 
ini,j (there is an algorithm A : ^ Q such that \A{i,j,n) — d{si,Sj)\ < 
2"";. 

5 is a numbered set, and the information that can be recovered from the 
numbers of ideal points is their mutual distances. Without loss of generality, we 
will suppose the numbering of S to be injective: it can always be made injective 
in an effective way. 
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Definition 5 We say that in a metric space {X, d), a sequence of points (x„)„gN 
converges recursively to a point x if there is an algorithm 13 : Q — N such that 
d{xn, x) < e for all n > Die). 

Definition 6 A point x & X is said to be computable if there is an algorithm 
.4 : N — > such that (^(?i))„gN converges recursively to x. 

We define the set of ideal balls to be B := {B{si, qj) : s,; e 5, < qj E Q} 
where B{x,r) ~ {y E X : d{x,y) < r} is an open ball. We fix a numbering 
B = {Bq, Bi, . . .} which makes the number of a ball effectively computable from 
its center and radius and vice versa. S is a countable basis of the topology. 

Definition 7 (Effective open sets) We say that an open set U is effective if 
there is an algorithm : N — > ,B such that U — IJ^^ A{n) . 

Observe that an algorithm which diverges on each input n enumerates the 
empty set, which is then an effective open set. Sequences of uniformly effective 
open sets are naturally defined. Moreover, if (Ui)i^fq is a sequence of uniformly 
effective open sets, then IJ- Ui is an effective open set. 

Definition 8 (Effective Gg-set) An effective G^-set is an intersection of a 
sequence of uniformly effective open sets. 

Obviously, an uniform intersection of effective G^-sets is also an effective 
G^-set. 

Let {X,Sx — {si ,S2 , ■■■},dx) and {Y,Sy = {sf , s|', ...}, rfy) be com- 
putable metric spaces. Let also Bf and BY be enumerations of the ideal balls 
in X and Y. A computable function X ^ Y is a function whose behavior can 
be computed by an algorithm up to any precision. For this it is sufficient that 
the pre-image of each ideal ball can be effectively enumerated by an algorithm. 

Definition 9 (Computable Functions) A function T : X ^ Y is com- 
putable if T^^{BY) is an effective open set, uniformly in i. That is, there is 
an algorithm A : N x N ^ B^ such that T^^{BY) = Un -^(^ Z^*" '^^^ 
A function T : X ^ Y is computable on D Q X if there are uniformly 
effective open sets such that T-'^{BJ) nD^UiOD. 

Remark 10 Intuitively, a function T is computable (on some domain C) if 
there is a computer program which computes T{x) (for x € C) in the following 
sense: on input e > 0, the program, along its run, asks the user for approxi- 
mations of X, and eventually halts and outputs an ideal point s € Y satisfying 
d{T(x),s) < e. This idea can be formalized, using for example the notion of 
oracle computation. The resulting notion coincides with the one given in the 
previous definitions. 

Recursive compactness is an assumption which will be needed in the follow- 
ing. Roughly, a compact set is recursively compact if the fact that it is covered 
by a finite collection of ideal balls can be tested algorithmically (for equivalence 
with the e-net approach and other properties of recursively compact set see [26]). 
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Definition 11 A set K C X is recursively compact if it is compact and 
there is a recursive function (/3 : N* — > N such that <yj(ii, . . . ,ip) halts if and only 
if {Bi^ ,Bi^) is a covering of K. 

3 Computing invariant measures 

3.1 Invariant measure and statistical properties 

Let X be a metric space, T : X ^-^ X a, Borel measurable map and fi a T- 
invariant Borel probability measure. A set A is called T-invariant if T~^{A) = 
A {mod 0). The system {X,T,fi) is said to be ergodic if each T-invariant set 
has total or null measure. In such systems the famous BirkhofF ergodic theorem 
says that time averages computed along /i-typical orbits coincides with space 
average with respect to /i. More precisely, for any / e L^{X, fi) it holds 



for ^ almost each x, where Sl^f + foT+... + fo T"-i. 

This shows that in an ergodic system, the statistical behavior of observables, 
under typical realizations of the system is given by the average of the observable 
made with the invariant measure. 

We say that a point x belongs to the basin of an invariant measure fi if (1) 
holds at X for each bounded continuous /. In case X is a, manifold (possibly with 
boundary) , a physical measure is an invariant measure whose basin has positive 
Lebesgue measure (for more details and a general survey see [47]). Computation 
of such measures will be the main subject of this section. 

3.1.1 The transfer operator 

A function T between metric spaces naturally induces a function Lt between 
probability measure spaces. This function Lt is linear and is called transfer 
operator (associated to T). Measures which are invariant for T are fixed points 
of Lt- 

Let us consider a computable metric space X and a measurable function 
T : X ^ X. Let us also consider the space PM{X) of Borel probability 
measures on X. 

Let us define the linear function Lt ■ PM{X) — > PM(X) by duality in the 
following way: if /i G PM{X) then LTijj) is such that 



The computation of invariant measures (and many other dynamical quan- 
tities) very often is done by computing the fixed points (or the spectrum) of 
this operator in a suitable function space. The most applied and studied strat- 
egy is to find a suitable finite dimensional approximation of Lt (restricted to 




(1) 
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a suitable function space) so reducing the problem to the computation of the 
corresponding relevant eigenvectors of a finite matrix. 

An example of this is done by discretizing the space X by a partition Ai and 
replacing the system by a (finite state) Markov Chain with transition probabil- 
ities 

~ m(A,) 

where m is the Lebesgue measure on X (see e.g. [20] [21] [37] ), then, taking 
finer and finer partitions it is possible to obtain in some cases that the finite 
dimensional model will converge to the real one (and its natural invariant mea- 
sure to the physical measure of the original system). In some case there is an 
estimation for this speed of convergence (see eg. [21] for a discussion), but a 
rigorous bound on the error (and then a real rigorous computation) is known 
only in a few cases (piecewise expanding or expanding maps, see [37]). 

Similar approaches consists in applying a kind of Faedo-Galerkin approxi- 
mation to the transfer operator by considering a complete Hilbert base of the 
function space and truncating the operator to the action on the first elements 
(see [46] ). 

Another approach is to consider a perturbation of the system by a small 
noise. The resulting transfer operator has a kernel. This operator then can be 
approximated by a finite dimensional one, again by Faedo-Galerkin method 
and relevant eigenvectors are calculated (see e.g. [17] [16]) then, if we prove that 
the physical measure of the original system can be obtained as a limit when 
the size of the noise tends to zero (this happen on uniformly hyperbolic system 
for example) we have a method which in principle can rigorously compute this 
measure. 

Variations on the method of partitions are given in [18, 19], while in [40] 
a different method, fastly converging, based on periodic points is exploited for 
piecewise analytic Markov maps. Another strategy to face the problem of com- 
putation of invariant measures consist in following the way the measure fi can 
be constructed and check that each step can be realized in an effective way. In 
some interesting examples we can obtain the physical measure as limit of iter- 
ates of the Lebesgue measure fi = lim„_j.oo U^{m) (recall that m is the Lebesgue 
measure). To prove computability of the main point is to explicitly estimate 
the speed of convergence to the limit. This sometimes can be done using the 
spectral properties of the system ([24]). 

Concluding, if the goal is to rigorously compute an invariant measure, most 
of the results which are in the today literature are partial. Indeed, beside 
being applied to a quite restricted class of systems, to compute the measure 
those methods need additional information. For example, the way to compute 
rigorously the finite dimensional approximation is often not done effectively, or 
the rate of convergence of the approximation is computed up to some constants 
depending on the system, which are not estimated. 

In the remaining part of the section we present some results, mainly from 
[26] explaining some result about rigorous computations of invariant measures. 
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These results have the advantage to give in principle an effective method for the 
rigorous computation of an invariant measure and give a quite general result, 
where all the needed assumptions are explicited. 

They have the disadvantage to not optimize computation time. So it is not 
clear if they can be implemented and used in practice. 

The rigorous framework into they are proved, however allows to see them 
as an investigation about the theoretical limits of rigorous computation of 
invariant measures. 

Moreover by this, also negative results can be proved. In particular, we give 
examples of computable systems having no computable invariant measures. 



3.2 Computability of measures 

In this section we explain precisely what we mean by computing a measure. 
This means, having an algorithm who is able to approximate the measure by 
"simple measures" up to any given approximation. 

Let us consider the space PAI{X) of Borel probability measures over X. Let 
Co{X) be the set of real- valued bounded continuous functions on X. We recall 
the notion of weak convergence of measures: 

Definition 12 is said to be weakly convergent to ii if J f d/i„ — > // dfi 
for each f € CqIx). 

Let us introduce the Wasserstein-Kantorovich distance between measures. 
Let and be two probability measures on X and consider: 



/ d/ii - / / d^ 



W^i(Mi,M2)= sup 

/ei-Lip(x) 

where l-Lip(X) is the space of functions on X having Lipschitz constant less 
than one. The distance Wi has the following useful properties 

Proposition 13 (see [1] Prop 7.1.5) 

1. Wi is a distance and if X is bounded, separable and complete, then PM{X) 
with this distance is a separable and complete metric space. 

2. If X is bounded, a sequence is convergent for the Wi metrics if and only 
if it is convergent for the weak topology. 

3. If X is compact PM(X) is compact with this metrics. 

Item (1) has an effective version: PM(X) inherits the computable metric 
structure of X. Indeed, given the set of ideal points of X we can naturally 
define a set of ideal points §pm(x) in PM{X) by considering finite rational 
convex combinations of the Dirac measures Sg supported on ideal points s € Sx ■ 
This is a dense subset of PM{X). The proof of the following proposition can 
be found in ([31]). 
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Proposition 14 If X bounded then {PM{X),Wi,^jPm{x)) is a computable 
metric space. 

A measure ^ is then computable if there is a sequence /x„ g ^pm(x) converg- 
ing recursively fast to in the Wi metric (and hence for the weak convergence). 

3.3 Computable invariant "regular" measures 

Invariant measures can be found as fixed points of the transfer operator, i.e. 
as solutions of the equation Wi{^,Lt{^j)) — 0. There is an abstract theorem 
allowing the computation of isolated fixed points (see [26] Corollary 3) of com- 
putable maps. To apply it is necessary to consider a space where the transfer 
transfer operator is computable (this space hopefully will contain the physical 
invariant measure). 

We remark that if T is not continuous then Lt is not necessarily continu- 
ous (this can be realized by applying Lt to some delta measure placed near a 
discontinuity point) hence not computable. Still, we have that Lt is continu- 
ous (and its modulus of continuity is computable) at all measures ^ which are 
"far enough" from the discontinuity set D. This is technically expressed by the 
condition ii{D) — 0. 

Proposition 15 Let X be a computable metric space and T : X X be a 
Junction which is computable on X \ D. Then Lt is computable on the set of 
measures 

PMd{X) := {/i e PM{X) : y.{D) = 0}. (2) 

From a practical point of view this proposition provides sufficient conditions 
to rigorously approximate the transfer operator by an algorithm. 

The above tools allow us to ensure the computability of Lt on a large class 
of measures. This will enable us to apply a general result on computation of 
fixed points and obtain 

Theorem 16 ([26], Theorem 3.2) Let X he a computable metric space and 
T a function which is computable on X \ D. Suppose there is a recursively 
compact set of probability measures V C PM{X) such that for every ^ ^ V , 
fJ.{D) — holds. Then every invariant measure isolated in V is computable. 
Moreover the theorem is uniform: there is an algorithm which takes as inputs 
finite descriptions ofT, V and an ideal ball in M(X) which isolates'^ an invariant 
measure /i, and outputs a finite description of fi. 

A trivial and general consequence of Theorem 16 is the following: 

Corollary 17 If a system as above is uniquely ergodic and its invariant measure 
jjL satisfies fJ-{D) = 0, then it is a computable measure. 

^If the invariant measure is unique in V the isolating ball is not necessary. 
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The main problem in the apphcation of Theorem 16 is the requirement that 
the invariant measure we are trying to compute be isolated in V. In general 
the space of invariant measures in a given dynamical system could be very large 
(an infinite dimensional convex subset of PM{X)) and physical measures have 
often some kind of particular regularity. To isolate a particular measure we can 
restrict and consider a subclass of "regular" measures. 

Let us consider the following seminorm: 

II II KB{x,r)) 

\\ma= sup . 

xex,r>o r"' 

If a and K are computable and X is recursively compact then 

Vo.,K^{^^ePM{x)■.\\^^\\^<K} (3) 

is recursively compact ([26]). This implies 

Proposition 18 Let X be recursively compact and T be computable on X\D , 
with dimniD) < oo. Then any invariant measure isolated in Va^K with a > 
dimniD) is computable. Once again, this is uniform in T,a,K. 

The above general propositions allow us to obtain as a corollary the com- 
putability of many absolutely continuous invariant measures. For the sake of 
simplicity, let us consider maps on the interval. 

Proposition 19 If X ^ [0, 1], T is computable on X \ D, with dimniD) < 1 
and {X, T) has a unique absolutely continuous invariant measure fi with bounded 
density, then fi is computable ( starting from T and a bound for the L°° norm 
of the invariant density). 

Similar results hold for maps on manifolds (again see [26]). 



3.3.1 Computing the measure from a description of the system in 
the class of piecewise expanding maps 

As it is well known, interesting examples of systems having a unique absolutely 
continuous invariant measure (with bounded density as required) are topologi- 
cally transitive piecewise expanding maps on the interval or expanding maps on 
manifolds. 

We show how to find a bound for the invariant density on piecewise expand- 
ing maps. By this, the invariant measure can be calculated starting from a 
description of T. 

Definition 20 A nonsingular function T : ([0, l],m.) — > ([0, l],77i) is said piece- 
wise expanding if 

•^For the seek of simplicity we will consider the simplest setting in which we can work and 
give precise estimations. Such class was generalized in several ways, we hence warn the rader 
that now in the current literature by piecewise expanding maps it it meant something slightly 
more general than the ones defined here. 
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1. There is a finite set of points di = 0, ^2, dn — 1 such that T\(^^. ,i.^^-j is 

and can he extended to a C'^map on [d^, d^+i] . 

2. inf(T') > 1 on the set where it is defined. 

3. T is topologically mixing. 

It is now well known that such maps have an unique ergodic invariant mea- 
sure with bounded variation density. 

Such density is also the unique fixed point of the (Perron Frobenius ) oper- 
ator^ L : Li[0, 1] ^ Li[0, 1] defined by 

We now show how to find a bound for such density, starting from the de- 
scription of the system, and then compute the associated invariant measure. 

The following proposition was proved in the celebrated paper [35] (Thm. 1 
and its proof) and it is now called Lasota-Yorke inequality. We give a precise 
statement where the constants are explicited. 

Proposition 21 LetT a piecewise hyperbolic map. If f is of bounded variation 
in [0, 1]. Let di, ...,dn be the discontinuity points ofT. 
IfX^ inf T'(x). Then 

a;G[04]-{(ii,...,d„} 

Var{Lf) < 2XVar{f) + B 



whe 



sup {\ (T'(i% \) 

^ xe[o,i]-{di,...,d„} ^ ^ " 
B = 



inf irTr)) min((ij - dj+i) ■ 

The following is an elementar fact about the behavior of real sequences 

Lemma 22 // a real sequence a„ is such that a„+i < Za„ -|- k for some I < 
1, fc > 0, then 

k 

sup(a„) < max(ao, j^^)- 

Proposition 23 // / is the density of the physical measure of a piecewise 
expanding map T as above and A > 2. Then 

Varif) < ^ 

where B is defined as above. 



*Note that this operator corresponds to the above cited transfes operator, but it acts on 
densities instead of measures. 
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Proof, (sketch) The topological mixing assumption implies that the map has 
only one invariant physical measures (see [45]). Let us use the above results 
iterating the constant density corresponding to the Lebesgue measure. Propo- 
sition 21 and Lemma 22 give that the variation of the iterates is bounded by 
jSix- Suppose that the limit measure has density /. By compactness of BV in 
L^, the above properties give a bound on the variation of / (see [35] proof of 
Thm. 1). ■ 

The following is a trivial consequence of the fact that ||/||oo < Var(f) + 
Jfdf,. 

Corollary 24 In the above situation ||/||oo < T~2x + ^■ 

The bound on the density of the invariant measure, together with Corollary 
19 gives the following uniform result on the computation of invariant measures 
of such maps. 

Theorem 25 Suppose a piecewise expanding map T and also its derivative is 
computable on [0, 1] — {di, ...,dn}, and also its extensions to the closed intervals 
[di,di+i] are computable. Then, the physical measure can be computed start- 
ing from a description of the system (the program computing the map and its 
derivative) . 

Proof, (sketch) Since T and T' are computable on each interval [di,di+i] we 
can compute a number A such that 1 < A < inf(r') (see [26] Proposition 3) If 
we consider the iterate T^instead of T the associated invariant density will be 
the same as the one of T, and if A^ > 2, then will satisfy all the assumptions 
needed on Proposition 23. Then we have a bound on the density and we can 
apply Corollary 19 to compute the invariant measure. ■ 



3.4 Unbounded densities, non uniformly hyperbolic maps 

The above results ensure computability of some absolutely continuous invariant 
measure with bounded density. If we are interested in situations where the 
density is unbounded, we can consider a new norm, "killing" singularities. 
Let us hence consider a computable function / : X — > M and 

II II f{x)fi{B{x,r)) 
y\fa= sup . 

If a and K are computable and X is recursively compact then 

V^.K = {tiePM{X):M^^^<K} (4) 

is recursively compact, and 18 also hold for this norm. If / is such that f{x) = 
when Muir-^Q — oo this can let the norm be finite when the density 

diverges. 
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As an example, where this can be apphed, let us consider the Manneville 
Pomeau type maps on the unit interval. These are maps of the type x — ^ 
X + x^(mod 1). 

When 1 < z < 2 the dynamics has a unique absolutely continuous invariant 
measure /x^. This measure has a density ez{x) which diverges at the origin, and 
ez{x) X x~^^^ and is bounded elsewhere (see [32] Section 10 and [45] Section 
3 e.g.). If we consider the norm with f{x) = x^ we have that ll/x^llj;^ is 

finite for each such z. By this it follows that for each such z the measure fi^ is 
computable. 

3.5 Computable systems without computable invariant mea- 
sures 

We have seen several techniques to establish the computability of many physical 
invariant measures. This raises naturally the following question: a computable 
systems does necessarily have a computable invariant measure? what about 
ergodic physical measures? 

The following is an easy example showing that this is not true in general 
even in quite simple systems, hence the whole question of computing invariant 
measures has some subtlety. 

Let us consider a system on the unit interval given as follows. Let r G (0, 1) 
be a lower semi-computable real number which is not computable. There is a 
computable sequence of rational numbers such that supj = r. For each i 
consider Ti{x) = max(a;,Ti) and 

T{x) ^Y.^-'T,. 
i>i 

The functions Ti are uniformly computable so T is also computable. 




'"I n ^3' ■ - r 



Figure 1: The map T. 
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The system ([0,1], T) is hence a computable dynamical system. T is non- 
decreasing, and T[x) > a; if and only if x < r. 

This system has a physical ergodic invariant measure which is 6t , the Dirac 
measure placed on r. The measure is physical because r attracts all the interval 
at its left. Since r is not computable then 5r is not computable. We remark 
that coherently with the previous theorems 5t is not isolated. 

It is easy to prove, by a simple dichotomy argument, that a computable 
function from [0, 1] to itself must have a computable fixed point. Hence it is not 
possible to construct a system over the interval having no computable invariant 
measure (we always have the 8 over the fixed point). With some more work 
we will see that such an example can be constructed on the circle. Indeed the 
following can be established 

Proposition 26 There is a computable, continuous map T on the circle having 
no computable invariant probability measure. 

For the description of the system and for applications to reverse mathematics 
see [26] (Proposition 12). 

4 Computability in Complex Dynamics 

Let X be a compact subset of the plane. Informally speaking, in order to draw 
the set K on a computer screen we need a program which is able to decide, 
given some precision e, if pixel p has to be colored or not. By representing pixel 
p by a ball B{p,e) where (p, e) e Q^, the question one would like to answer is: 
does B{p,e) intersects K ? The following definition captures exactly this idea: 

Definition 27 A compact set K cM.'^ is said to be computable if there is an 
algorithm A such that, upon input {p, e) .■ 

• halts and outputs "yes" if K O B{p, e) ^ 0, 

• halts and outputs "no" if K D B{p,e) = 0, 

• run forever otherwise. 

We remark that if a compact set is computable under the definition above, 
then it is recursively compact in the sense of Definition 11. The converse is 
however false. 

4.1 Julia sets and the Brolin-Lyubich measure 

The question of whether Julia sets are computable under this definition has been 
completely solved in a series of papers by Binder, Braverman and Yampolsky. 
See [4, 5, 11, 12, 13]. Here we review some of their results. For simplicity, we 
give the definition of the Julia set only for quadratic polynomials. Consider a 
quadratic polynomial 

Pc{z) = + c : C ^ C. 
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Obviously, there exists a number M such that if |2;| > M, then the iterates 
of z under P will uniformly scape to oo. The filled Julia set is then defined 
by: 



i^c = {2 e C sup |P"(z)| < 00}. 



That is, the set of points whose orbit remains bounded. For the filled Julia 
set one has the following result: 

Theorem 28 ([13]) The filled Julia set is always computable. 
The Julia set can be now defined by: 



where d{A) denotes the boundary of A. The Julia set is the repeller of the 
dynamical system generated by Pc- For all but finitely many points, the limit of 
the n-th preimages P~"(z) coincides with the Julia set Jc. The dynamics of Pc 
on the set Jc is chaotic, again rendering numerical simulation of individual orbits 
impractical. Yet Julia sets are among the most drawn mathematical objects, 
and countless programs have been written for visualizing them. 
In spite of this, the following result was shown in [11]. 

Theorem 29 There exist computable quadratic polynomials Pc{z) — z'^ + c 
such that the Julia set Jc is not computable. 

This phenomenon of non-computability is rather subtle and rare. For a 
detailed exposition, the reader is referred to the monograph [12]. 

Thus, we cannot accurately simulate the set of limit points of the preimages 
{Pc)~'^{z), but what about their statistical distribution? The question makes 
sense, as for all z 7^ 00 and every continuous test function ip, the averages 



where A is the Brolin-Lyubich probability measure [14, 39] supported on the Julia 
set Jc- We can thus ask whether the the Brolin-Lyubich measure is computable. 
Even if Jc = Supp(A) is not a computable set, the answer does not a priori have 
to be negative. In fact, the following result holds: 

Theorem 30 ([6]) The Brolin-Lyubich measure is always computable. 

The proof of the previous result does not involve (perhaps surprisingly) much 
analytic machinery, but it follows from some general principles in the same spirit 
of the ideas introduced in Section 3.3. 



Jc = dK, 




19 



4.2 The Mandelbrot set 



Another important object in complex dynamics is the Mandelbrot set, which is 
defined as follows: 

M:={ceC:Oe K^}. 

That is, the set of parameters of quadratic polynomials for which the orbit 
of remains bounded. As for Julia sets, there exists many computer programs 
to visualize it, and the question arise whether this set is actually computable 
in a rigorous sense. This is, up to date, an open question. However, some 
partial results have been obtained. For instance, in [29] it is proved that the 
Mandelbrot set is recursively compact in the sense of Definition 11. There exists 
also analytical questions about the Mandelbrot set that remain unsolved. For 
instance, it is unknown whether it is locally connected. It has been conjectured 
that this is indeed the case. Computability of the Mandelbrot set is closely 
related to this conjecture, as it has been observed in [29]: 

Theorem 31 ([29]) // the Mandelbrot set is locally connected, then it is com- 
putable. 

It is known that local connectivity implies another well known conjecture, 
namely the hyperbolicity conjecture, which in turn also implies computability. 
For a detailed exposition in this subject the reader is referred to [15]. 

5 Computing the speed of ergodic convergence 

As recalled before, the Birkhoff ergodic theorem tells that, if the system is 
ergodic, there is a full measure set of points for which the averages of the values 
of the observable / along its trajectory (time averages) coincides with the spatial 
average of the observable /. Similar results can be obtained for the convergence 
in the norm, and others. Many, more refined results are linked to the speed 
of convergence of this limit. And the question naturally arise, if there is a 
possibility to compute this speed of convergence in a sense similar to Definition 
5. 

In the paper [2] some abstract results imply that in a computable ergodic 
dynamical system, the speed of convergence of such averages can be algorithmi- 
cally estimated. On the other hand it is also shown that there are non ergodic 
systems where this kind of estimations are not possible. In [27] a very short 
proof of this result for ergodic systems is shown. We expose the precise result 
(Theorem 39) and the short proof in this section. 

5.0.1 Convergence of random variables 

We first precise what is meant by "compute the speed" in some pointwise a.e. 
convergence. 
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Definition 32 A random variable on {X, /i) is a measurable function f : 
X ^R. 

Definition 33 Random variables fn effectively converge in probability to 

f if for each e > 0, fi{x : \ fn{x) — f{x)\ < e} converges effectively to I, uniformly 
in e. That is, there is a computable function n(e, S) such that for all n > n(e, S), 
f^{\fn-f\>e}<S. 

Definition 34 Random variables fn effectively converge almost surely to 

f */ fn = supj.>„ l/fc — /I effectively converge in probability to 0. 

Definition 35 A computable probability space is a pair {X, /i) where X 
is a computable metric space and /i a computable Borel probability measure on 
X. 

Let y be a computable metric space. A function {X, /i) -> F is almost 
everywhere computable (a.e. computable for short) if it is computable on 
an effective G^-set of measure one, denoted by dom/ and called the domain of 
computability of f. 

A morphism of computable probability spaces / : {X, n) — )■ (Y, v) is a 
morphism of probability spaces which is a.e. computable. 

Remark 36 A sequence of functions /„ is uniformly a. e. computable if the 
functions are uniformly computable on their respective domains, which are uni- 
formly effective Gg-sets. Observe that in this case, intersecting all the domains 
provides an effective Gg-set on which all fn are computable. In the following we 
will apply this principle to the iterates fn — of an a. e. computable function 
T : X ^ X, which are uniformly a.e. computable. 

Remark 37 The space L^{X,ii) (resp. L'^(X,ii)) can be made a computable 
metric space, choosing some dense set of bounded computable functions as ideal 
elements. We say that an integrable function f : X ^ M. is L^{X, ^) -computable 
if its equivalence class is a computable element of the computable metric space 
L^(X, /i). Of course, if f ^ g ^i-a.e., then f is L^{X, computable if and 
only if g is. Basic operations on L^{X,^), such as addition, multiplication by 
a scalar, min, max etc. are computable. Moreover, if T : X ^ X preserves jjl 
and T is a.e. computable, then / — > / o T (from to ) is computable (see 
[30]). 

Let us call {X, /i, T) a computable ergodic system if {X, /i) is a com- 
putable probability space where T is a measure preserving morphism and {X, fi, T) 
is ergodic. Let ||/|| denote the norm or the norm. 

Proposition 38 Let {X, ii, T) be a computable ergodic system. Let f be a com- 
putable element of L^{X,^) (resp. L'^{X,fj,)). 

The L^ convergence ( resp. L^ convergence ) of the Birkhoff averages An = 
{f+foT+. . .'^foT"'~^)/n is effective. That is: there is an algorithmn : Q — )■ N 
such that for each m > n(e) ||A„j — J /d/i|| < e. Moreover the algorithm depends 
effectively on T, ^, f. 
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Proof. Replacing / with f — J fdn, we can assume that J fdfi — 0. The 
sequence ||A„|| is computable ( see Remark 37 ) and converges to by the 
ergodic theorems. 

Given p d N, we write m d N as m = rip + k with < k < p. Then 



Anp+k = -^{TpApoTP' + kAkoTP- 
np + k ^ ^-^ 



.1=0 



n 

< M + M. 

Let e > 0. We can compute some p = p{e) such that \\Ap\\ < e/2. Then we 
can compute some n(e) > |||/||. The function m(e) :— n{e)p{e) is computable 
and for all m > m(e), \\Ajn\\ < e. D 

5.1 Effective almost sure convergence 

Now we use the above result to find a computable estimation for the a.s. speed 
of convergence. 

Theorem 39 Let (X, /i,T) be a computable ergodic system. If f is L^{X,^)- 
computable, then the a.s. convergence is effective. Moreover, the rate of con- 
vergence can be computed as above starting from T,fi, f. 

This will be proved by the following 

Proposition 40 // / is L^{X, ji)- computable as above, and ||/||oo is bounded, 
then the almost-sure convergence is effective (uniformly in f and a bound on 
and on T , /i ). 



To prove this we will use the Maximal ergodic theorem: 
Lemma 41 (Maximal ergodic theorem) For f E L^{X,ii) and S > 0, 

f^{{sup\Ai\>S})<-\\f\U. 

n 

The idea is simple: compute some p such that | |v4^| |i is small, apply the max- 
imal ergodic theorem to g := A^, and then there is no, that can be computed, 
such that A(^ is close to A^^ for n> uq. 

Proof (of Proposition 40) Let e,S > 0. Compute p such that ||^^|| < Se/2. 
Applying the maximal ergodic theorem to g := gives: 

^,i{snp\Ai\>6/2})<e. (5) 



22 



Now, A-l^ is not far from A(^: expanding A-^^, one can check that 

where u = (p - l)f + {p~ 2)f oT + . . . + f oTP-\ \\u\\^ < ^^||/|U so if 
n>no> 4(p- l)||/||oo/'5, then \\Ai - AfJ^ < S/2. As a resuh, if \Alix)\ > S 
for some n > uq, then |Afj(a::)| > S/2. From (5), we then derive 

M({sup \Ai\>S})<e. 

n>no 

As uq can be computed from S and e, we get the result. □ 

Remark 42 This result applies uniformly to a uniform sequence of computable 
L°°{X,^) observables fn- 

We now extend this to L^{X, /i)-computable functions, using the density of 
L°°(X,y) in L^{X,ii). 

Proof, (of Theorem 39) Let e,(5 > 0. For M e N, let us consider f'f^,j G 
L°° (X, pi) defined as 



min(/,M) if f{x)>0 
max(/,-Af) iff{x)<0. 



Compute M such that ||/— /mHi < 5e. Applying Proposition 40 to /]^^ 
gives some no such that Ai({sup„>„„ \ An"\ > 5}) < e. Applying Lemma 41 to 

fli^f- I'm gives M({sup„ > S}) < e. As a result, /i({sup„>„^ > 

2d}) < 2e. a 

Remark 43 We remark that a bounded a.e. computable function, as defined 
in Definition 35 is a computable element of L^{X,fj,) (see [30]). Conversely, 
if f is a computable element of L^{X,^) then there is a sequence of uniformly 
computable functions fn that effectively converge pL-a.e. to f . 



6 Computing pseudorandom points, construc- 
tive ergodic theorem 

Let X again be a computable metric space and a computable probability 
measure on it. Suppose X is complete. 

Points satisfying the above recalled pointwise ergodic theorem for an observ- 
able /, will be called typical for f . 

Points which are typical for each / which is continuous with compact support 
are called typical for the measure /i ( and for the dynamics) . 

The set of computable points in X (see Definition 6) being countable, is 
a very small (invariant) set, compared to the whole space. For this reason, a 
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computable point could be rarely be expected to be typical for the dynamics, 
as defined before. More precisely, the Birkhoff ergodic theorem and other the- 
orems which hold for a full measure set, cannot help to decide if there exist a 
computable point which is typical for the dynamics. Nevertheless computable 
points are the only points we can use when we perform a simulation or some 
explicit computation on a computer. 

A number of theoretical questions arise naturally from all these facts. Due to 
the importance of forecasting and simulation of a dynamical system's behavior, 
these questions also have some practical motivation. 

Problem 44 Since simulations can only start with computable initial condi- 
tions, given some typical statistical behavior of a dynamical system, is there 
some computable initial condition realizing this behavior? how to choose such 
points? 

Such points could be called pseudorandom points, and a result showing its 
existence and their computability from the description of the system could be 
seen as a constructive version of the pointwise ergodic theorem. 

Meaningful simulations, showing typical behaviors of the dynamics can be 
performed if computable, pseudorandom initial conditions exist ( and can be 
computed from the description of the system). Thanks to a kind of effective 
Borcl-Cantclli lemma, in [24] the above problem was solved affirmatively for 
a class of systems which satisfies certain technical assumptions which includes 
systems whose decay of correlation is faster that C\og^{time). After the results 
on the estimation of the rate of convergence given in the previous section we can 
remove the technical assumption on the speed of decay of correlations. We will 
moreover show how the result is uniform also in T and /i (the pseudorandom 
points will be calculated from a description of the system). 

The following result ([24], Theorem 2 or [27]) shows that given a sequence 
/„ which converges effectively a.s. to / and given its speed of convergence then 
there it is possible to compute points Xi for which fn{xi) — f{xi). 

Theorem 45 Let X be a complete metric space. Let /«,/ be uniformly a.e. 
computable random variables. If /„ effectively converges almost surely to f then 
the set {x : fn{x) — > f{x)} contains a sequence of uniformly computable points 
which is dense in Supp(/i). Moreover, the effective sequence found above depends 
algorithmically on fn and on the function n{5, e) giving the rate of convergence 

Since by the results of the previous section, n{S, e) can be calculated starting 
from r, ^ and /. This result hence can be directly used to find typical points 
for the dynamics. Indeed the following holds (see [27] for the details) 

Theorem 46 // (X, /i,T) is a computable ergodic system, f is L^(X,fj,) and 
a. e. computable then there is a uniform sequence Xi of computable points which 
is dense on the support of /i such that for each i 
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Moreover this sequence can be computed starting from a description oJT, /i and 
/• 

Since the above result is uniform in / and it is possible to construct a r.e. set 
of computable observables which is dense in the space of compactly supported 
continuous functions we can also obtain the following (see again [27] for the 
details) 

Theorem 47 // (X,fi,T) is a computable ergodic system then there is a uni- 
form sequence x„ of computable points which is dense on the support of fi such 
that for each n, Xn is fi— typical. Moreover this sequence can be computed start- 
ing from a description of Tand fj. 

Now, as an application we get together the uniform results about computa- 
tion of invariant measures for the class of piecewise expanding functions shown 
in Subsection 3.3.1 to show that in that class the pseudorandom points can be 
calculated starting from a description of the map. Indeed since in that class the 
physical invariant measure /i can be calculated starting from T (see Theorem 
25). We hence obtain 

Corollary 48 Each piecevise expanding map, with computable derivative, as in 
the assumptions of Theorem 25 has a sequence of pseudorandom points which 
is dense in the support of the measure. Moreover this sequence can be computed 
starting from the description of the system. 

Hence establishing some kind of constructive versions of the ergodic theorem. 

7 Conclusions and directions 

In this article we have reviewed some results about the rigorous computation of 
invariant measures, invariant sets and typical points. Here, the sentence rigorous 
computation means "computable by a Turing Machine" . Thus, this can be seen 
as a theoretical study of which infinite objects in dynamics can be arbitrarily 
well approximated by a modern computer (in an absolute sense), and which 
cannot. In this line, we presented some general principles and techniques that 
allow the computation of the relevant objects in several different situations. On 
the other hand, we also presented some examples in which the computation of 
the relevant objects is not possible at all, stating some theoretical limits to the 
abilities of computers when used to simulate dynamical systems. 

The examples of the second kind, however, seem to be rather rare. An 
important question is therefore whether this phenomenon of non-computability 
is robust or prevalent in any sense, or if it is rather exceptional. For example, 
one could ask whether the non-computability is destroyed by small changes 
in the dynamics or whether the non-computability occurs with small or null 
probability. 

Besides, in this article we have not considered the efficiency (and therefore 
the feasibility) of any of the algorithms we have developed. An important (and 
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more difficult) remaining task is therefore the development of a resource bounded 
version of the study presented in this paper. In the case of Julia sets, for 
instance, it has been shown in [9, 43, 10] that hyperbolic Julia sets, as well as 
some Julia sets with parabolics, are computable in polynomial time. On the 
other hand, in [5] it was shown that there exists computable Siegel quadratic 
Julia sets whose time complexity is arbitrarily high. 

For the purpose to compute the invariant measure form the description of 
the system, in section 3.3.1 we had to give explicit estimations on the constants 
in the Lasota Yorke inequality. This step is important also when techiques 
different from ours are used (see [37] e.g.). Similar estimations could be done 
in other classes of systems, following the way the Lasota Yorke inequality is 
proved in each class (although, sometimes this is not a completely trivial task 
and requires the "effectivization" of some step in the proof). A more general 
method to have an estimation for the constants or other ways to get information 
on the regularity of the invariant measure would be useful. 
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